arXiv:l502.02774v2 [cond-mat.mes-hall] 21Jul2015 


Enhancement of tunneling density of states at a Y junction of 
spin-1 Tomonaga-Luttinger liquid wires 

Sougata Mardanya 1 and Amit Agarwal 1 ’* 

1 Department of Physics , Indian Institute of Technoloqy Kanpur, Kanpur - 208016, India 

(Dated: July 22, 2015) 

We calculate the tunneling density of states (TDOS) in a dissipationless three wire junction of 
interacting spin-1/2 electrons, and find an anomalous enhancement of the TDOS in the zero bias 
limit, even for repulsive interactions for several bosonic fixed points. This enhancement is physically 
related to the reflection of holes from the junction for incident electrons, and it occurs only in the 
vicinity of the junction (x < v m m/2w where v m i n is the minimum of the velocity of charge or spin 
excitations and u is the bias frequency), crossing over to the bulk value which is always suppressed, 
at larger distances. The TDOS exponent can be directly probed in a STM experiment by measuring 
the differential tunneling conductance as a function of either the bias voltage or temperature as 
done in C. Blumenstein et al, Nat. Phys. 7, 776 (2011). 
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I. INTRODUCTION 

One dimensional (ID) quantum wires with strong 
electron-electron (e-e) interactions are described by the 
Tomonaga-Luttinger liquid (TLL) theory [1-6], in which 
the low energy excitations are collective density oscil¬ 
lations. These density oscillations or plasmon modes, 
are fundamentally different from their 2D and 3D coun¬ 
terpart, i.e., Landau’s quasi-particle excitations, which 
are described very successfully by the Fermi liquid the¬ 
ory [7]. This leads to unique physics in ID, such as 
the spin-charge separation [8 and 9], the phenomena of 
charge fractionalization [10-15], power law behaviour of 
the differential tunneling conductance (with bias voltage 
or temperature) in quantum Hall edge states [16], carbon 
nanotubes [17], tailored quantum wires in semiconductor 
heterostructures [18] etc. TLL behaviour has also been 
observed in self-aligned Au atomic chains of single-atom 
width on germanium surface using scanning tunneling 
spectroscopy and photoemission [19]. 

In this paper we explore the local single particle tun¬ 
neling density of states (TDOS) at a junction of three 
spin full TLL wires, which can be experimentally probed 
in a scanning tunneling microscope (STM) experiment 
[20]. Such Y-junctions have already been realized exper¬ 
imentally in crossed single-walled carbon nanotubes [21 
and 22], and have been explored very actively in the lit¬ 
erature [23-40]. The differential tunneling conductance 
of a STM tip maintained at a finite bias voltage with 
respect to the wire [see Fig. 1] is a direct probe of the 
locally available electronic states in the TLL, assuming 
the density of states in the STM tip to be a constant. In 
a TLL wire the STM current has been shown to vary as 
a power of the bias voltage: dl/dV oc p(u) ex co A ~ 1 , with 
the TDOS exponent A depending on the strength of the 
e-e interaction strength [41-47]. The TDOS exponent in 
a spinless two wire junction tuned to the connected (or a 
single wire without impurity) [5 and 48] and disconnected 
fixed points (single wire with impurity) [5, 49-51], are 


known to be A = {g + g 1 )/2 and A = 1 /g respectively, 
where g denotes the TLL interaction parameter. 

Earlier theoretical studies of Y-junctions have used the 
fermionic language and the weak interaction renormaliza¬ 
tion group (RG) approach [24] , or the bosonic and confor¬ 
mal field theory language [26 and 28] , or other numerical 
methods such as the functional RG [27], and were pri¬ 
marily focused on the stability analysis of various fixed 
points of the junction and the associated DC conductiv¬ 
ity. A detailed analysis of the fixed points of a three 
wire junction formed from spinless interacting electrons 
was done in Ref. [26] which was later extended to include 
spin-1/2 electrons giving a much richer phase diagram in 
the parameter space of spin and charge interactions [31]. 

The tunneling density of states in a three wire Y- 
junction of spinless electrons was explored in Ref. [33], 
where the authors highlighted an anomalous enhance¬ 
ment of the TDOS in the vicinity of the junction for 
several bosonic fixed points. Here we consider a more re¬ 
alistic case of a three wire junction formed from spin-1/2 
interacting electrons, and obtain an analytic expression 
for the TDOS in the vicinity of the Y-junction. This pa¬ 
per is organized as follows: In Sec. II we describe the de¬ 
tails of the spin-1/2 TLL Y-junction and show that both 
the e-e interactions in the wire and the current partition¬ 
ing boundary conditions at the junction, can be treated 
using bosonization with delayed evaluation of the bound¬ 
ary conditions [26] taking advantage of the spin charge 
separation in ID. In Sec. Ill, we calculate the TDOS of 
the spin-1/2 TLL Y-junction which shows an anomalous 
enhancement for several fixed points even for repulsive 
interactions (see Figs. 2 and 3). Finally we summarize 
our findings in Sec. IV. 


II. MODEL AND BOSONIZATION 

In this section we review the TLL description of in¬ 
teracting ID wires with spin-1/2 electrons which gives 
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FIG. 1. Schematic of a Y junction of three spin-1/2 TLL 
wires, and a STM tip for measuring the TDOS. The red and 
green lines indicate the charge and spin sectors respectively, 
and the inter-wire tunneling region is represented by the cir¬ 
cular region. 


rise to spin charge separation in the bulk of the wires. 
The dissipation-less junction is then described by means 
of additional tunneling between the wires, which leads 
to spin and charge current partitioning boundary condi¬ 
tion at the junction. The system of three spin-1/2 TLL 
wires with the junction (Y-junction) is then explored us¬ 
ing the delayed evaluation of boundary conditions first 
developed in Ref. [26]. We note that the technique of 
delayed evaluation of boundary conditions is an exten¬ 
sion of the plasmon scattering approach first developed 
in Refs. [10-12], where the outgoing chiral currents are 
linearly related to the incoming currents and the conduc¬ 
tance is related to the scattering of plasmons. 


A. Bosonization of the wires 

Let us consider the Y-junction to be made of three 
identical, semi-infinite (x > 0) TLL wires with spin-| 
electrons, with the junction being located at the origin. 
For simplicity we consider all the wires to have the same 
short ranged e-e interaction strength, and same Fermi ve¬ 
locity for both up-spin and down-spin electrons in all the 
wires. In each wire, the physical fermionic field for each 
spin species a can be expressed in terms of slowly varying 
incoming (left movers) and outgoing (right movers) chi¬ 
ral fields as ip a (x) = e lkFX xl>° + e ~ lkFX %/>£. The fermionic 
field of an outgoing (right-movers) and incoming (left- 
movers) electron with spin a can be bosonized as, 

V>°(x) = —=F 0 a e i2wNo ^ x ~ vt) / L e ~ iv °( x) , 

V27ra 

%l>,i{x ) = ~^=F IrJ e O*I\Mx+i;t)/£ , (1) 

V27TOI 

where with a = O/I indicating the propaga¬ 

tion direction, represents the collective density excita¬ 
tions in ID. In Eq. (1) F aa denotes the Klein factor (or 
ladder operators) which increase and decrease the corre¬ 
sponding fermion number and satisfy anti-commutation 
relations among themselves, a is the inverse ultraviolet 


(short-distance) cut off and N aa counts the number of in¬ 
coming or outgoing chiral particles of spin a with respect 
to the filled Fermi sea. The bosonic field operators can be 
expressed in terms of bosonic creation and annihilation 
operators as, 



and the bosonic fields satisfy the following equal 
time commutation relations in a non-interacting theory, 
[<^“(x), </?“', (a/)] = ±?7r0(x — x')5 aa i5 aa ', where the + (-) 
sign arises for a denoting the outgoing (incoming) mode 
in each wire, and @(x) is the Heaviside step function. 

The corresponding incoming and outgoing charge den¬ 
sity and the current fields in each wire are given by 


o _ TO _ d tFa 

Pa 2tt ’ a 2tt 

/ _ d x (pi j _ d t <pj 

Pa 2tt ’ a 2tt 


( 3 ) 


Anticipating spin-charge separation (decoupling) in ID, 
we define the independent charge and spin fluctuation 
fields in each wire which commute with each other: 




( 4 ) 


Using the notation, v = c/s to represent either charge 
or spin degrees of freedom, these can in turn be used to 
form the pair of canonically conjugate phase fields in each 
wire: 


= 2= [( p° + ipl] , 9 V = 2= [<pl - <p°] , (5) 

which, in a non-interacting theory satisfy the commu¬ 
tation relation = — ittQ(x — x'). The 

Hamiltonian for each wire, including the effect of short 
ranged e-e interaction is now expressed in terms of these 
fields as 



g v (y<t> v ) 2 


-(V1U) 2 


( 6 ) 


where v c and v s denote the charge and spin velocities and 
g c and g s are the two interaction parameters for charge 
and spin sector, respectively, [6] which we assume to be 
the same in all the wires. Note that g c / s = 1 corresponds 
to the noninteracting case, 0 < g c < 1 to repulsive e-e in¬ 
teractions, and 5 S 1 to a broken spin SU(2) symmetry. 
In the absence of an external magnetic field, or spin- 
dependent interactions, g s = 1 due to the underlaying 
SU(2) symmetry in the spin space. In Eq. (6), L denotes 
the length of the wires, which is assumed to be much 
longer than the width of the low-energy wave packets, 
LkF 1, which allows us to use a continuum descrip¬ 
tion. Furthermore we have completely ignored phonons 
and disorder in the TLL wires. 
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The Hamiltonian in Eq. (6), which includes e-e interac¬ 
tions, can be diagonalized in terms of free bosonic fields in 
each wire, which satisfy the equal time canonical commu¬ 
tation relation: [</>“ (a;),/>“, (a/)] = ai ttQ(x — x')6 aa i5 aa i. 
This is achieved by scaling the interacting fields in Eq. (6) 
as 

<t> v = , and 0 V = g~ 0 U . (7) 

The incoming and outgoing interacting helds in each 
wire, can now be expressed in terms of the scaled 
free helds through the usual Bogoliubov transformation, 
given by 

<P% = [(! + 9cW c + (1 - <7c)<^c] 

+ a 9 V ^ = + 9s)<fis + (1 _ 3s)v?“] • (8) 

In the above equation a is dehned such that if a = I 
then a = O and vice versa and a = ±1 for up-spin and 
down-spin electrons respectively (when not used as a sub¬ 
script). 

B. Bosonization of the junction 

In addition to the Hamiltonian in Eq. (6) which de¬ 
scribes each of the three disconnected spin-1/2 inter¬ 
acting wires, to describe the junction we need to im¬ 
pose additional boundary conditions on the helds in each 
wire at the junction. As a consequence of the spin- 
charge separation in a spin-1/2 TLL wire, the charge 
and the spin sectors satisfy independent boundary condi¬ 
tions. Following the standard procedure [26, 28, 31, and 
33], we consider a current splitting matrix M, which re¬ 
lates the incoming charge and spin currents to the out¬ 
going spin and charge currents (in the noninteracting 
theory), and consequently the incoming and outgoing 
bosonic helds, i.e, ^ , which leads 

to <j>2 ti = 4>l j (where v = c/s denotes the 

charge or the spin sector). Here we have neglected an 
integration constant which plays no physical role in the 
computation of TDOS exponent or in the scaling dimen¬ 
sions of various operators. To enforce the condition that 
the boundary condition specihed by the matrix M repre¬ 
sents a hxed point of the theory, the incoming and out¬ 
going helds must satisfy the usual bosonic commutation 
relations. This restricts the matrix M to be orthogonal 
and the orthogonality condition of M also implies that 
there is no dissipation at the junction [34]. Furthermore, 
the conservation of charge and spin current at the junc¬ 
tion, ensure that each of the rows of the matrix M add 
up to unity. 

For a junction of two or more TLL spin-1/2 wires, one 
has to impose independent boundary conditions on the 
charge and spin sectors, and thus all hxed points will be 
represented by a product of two matrices: M C M S , with 
the hrst matrix specifying the boundary for the charge 


sector, and the second matrix specifying the boundary 
condition for the spin sector. For the case of a two-wire 
junction, there are only two possibilities for the M ma¬ 
trices which are given by, 

Rn =(o!) and Rd = (!!)- ( 0 ) 

with Rn representing the disconnected or reflecting hxed 
point (the subscript N denotes Neumann boundary con¬ 
dition) and Rjj representing the connected or transmit¬ 
ting hxed point respectively (the subscript D stands for 
Diriclrlet boundary condition), with both of them pre¬ 
serving time reversal symmetry (TRS). Now, all possible 
hxed point or boundary conditions for a two wire junction 
of spin-1/2 electrons are given by RnRn (both charge 
and spin sectors disconnected), RnRd (charge sector 
disconnected and spin sector connected), R^Rd (both 
charge and spin sectors connected) and finally RjjRn 
(charge sector connected and spin sector disconnected). 

For a three wire charge and spin conserving Y junc¬ 
tion, all current splitting orthogonal matrices M can be 
parametrized by a single continuous parameter 0 as in 
the case of spin-less electron wires [28 and 33] . For spin- 
1/2 wires, there are two such continuous parameters 0 C 
and 9 S , specifying the boundary condition (or hxed point) 
corresponding to the charge and spin degree of freedom 
respectively. The 3x3 matrices M specifying the bound¬ 
ary condition for each sector, preserve time reversal sym¬ 
metry (TRS), only if they are symmetric, and based on 
this they can be divided into two classes: det(Mi„) = 1 
and det(M 2 „) = — 1. These two classes have a form ex¬ 
plicitly given by 

( a b c\ (b a c\ 

cab J , Ma* = j a c b J , (10) 

b c a) v \c b a) v 

where a = (l + 2cos6 | 1/ )/3, b = (1 — cos0„ + -\/3sin0„)/3, 
and c = (1 — cos 6 V — v / 3sin0^)/3. This implies four 
families of hxed points for a Y-junction of interacting 
spin-1/2 wires: Mi c Mi s , Mi c M 2s , M 2 C Mi s and finally 
M 2 c M 2s . Of these only the M lc M ls corresponds to a 
symmetric (in the wire indices) class of hxed points, 
while the other three specify an asymmetric class of hxed 
points (with broken symmetry in the wire indices). 

The matrix M connects the incoming and outgoing 
held in a non-interacting theory. In presence of e-e inter¬ 
actions in the wire, it also has to undergo a Bogoliubov 
transformation: M -> M so that it connects the incom¬ 
ing and outgoing effective ‘free’ helds at the junction via 
the relation, </>[A(0,t) = JY(M„)ij <^(0, t). The Bo¬ 
goliubov transformed matrix is given by [28 and 33], 

M„ = [(1 + ^)1- (l-^)M 1/ ] _1 [(l + g 1 ,)M y - (1-^)1] . 

(11) 

The M 2l y class of matrices have an interesting property: 
(M 2l ,) 2 = I. Asa consequence M 2; , = M 2 „, which implies 
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that at the junction both interacting and free fields have 
identical properties. For Mij,, the matrix Mi„ still has 
the same form as Mi, but with the corresponding ma¬ 
trix elements given by a = (3 g 2 — 1 + (3 g 2 + 1) cos 9 v )/b 
and b/c = 2(1 — cos 9„ ± y/3g„ sin9 u )/5, where S = 
3[l + <^-|-(<^ — 1) cos$„]. Note that the formulation of de¬ 
layed evaluation of boundary condition described in this 
section and Eq. (11) is applicable to any number of inter¬ 
acting one-dimensional wires connected at a dissipation¬ 
less junction described by a bosonic fixed point. 

A detailed and systematic study of the stability of var¬ 
ious fixed points for a Y-junction of spin-1/2 electrons 
has been done in Ref. [31], using conformal field the¬ 
ory as well as bosonization with delayed evaluation of 
boundary conditions, in the g s — g c parameter space and 
several important fixed points with a unique attractive 
basin have been identified. For the sake of completeness, 
and as a check of our calculations, we report the scaling 
dimensions of all the spin-preserving single particle and 
pair tunneling operators, for all possible fixed points in 
Appendix. 

For the Mi c Mi s class of fixed points, 9 C / S = n 
or [a,b,c] = [—1/3,2/3,2/3], represents the so called 
DD (Dirichlet-Dirichlet) fixed point for charge/spin sec¬ 
tor, 9 C / S = 0 or [a, b, c) = [1,0,0] indicates the NN 
(Neumann-Neumann) or disconnected fixed point in 
which there is no tunneling between any pair of wires. 
The case of 9 C / S = ([ a,b,c ] = [0,1,0]) and 9 C / S = ^ 

([a, b, c] = [0,0,1]) corresponds to the chiral x _ X~ and 
X + X + fixed points, respectively. Other important fixed 
points belonging to the Mi c Mi s class with a unique at¬ 
tractive basin are specified as DN for 9 C = tt and 9 S = 0 
and ND for 9 C = 0 and 9 S = tt. 

For the M 2 C M 2 S class of fixed points, 9 C / S = 0, 27t/ 3 or 
4 -7 t/3 represents the asymmetric set of fixed points called 
DaDa , where two of the three wires are connected in 
both the charge and spin sectors while one of the wires 
is completely disconnected. All the three cases are iden¬ 
tical and we will discuss the case of 9 C = 9 S = 0, where 
wires 1 and 2 are connected in both the charge and spin 
sectors and wire 3 is completely disconnected. Note that 
even though the DaDa fixed point is asymmetric in the 
wire indices, it preserves TRS since only two of the wires 
are connected. There are several other interesting fixed 
points as well, however only the ones discussed here have 
a unique basin of attraction in the g s — g c plane. Along 
the SU(2) invariant line, g s = 1, the NN fixed point 
is stable for g c < 1 , the chiral fixed point is stable for 
1 < 9c < 3, and the DN stable fixed point is stable for 
5c > 3 [31]. 

Finally, we note that the matrices M„ specifying the 
boundary condition for the charge and spin degrees of 
freedom at the junction are also associated with the 
charge and spin conductance tensors associated with each 
fixed point. The outgoing charge and spin current in 
wire j is given by // = G'/.V/', where the superscript 
v = c/s, and V£ is the voltage applied on wire fc, and 
Vjf = fj,kt — Pkl is the chemical potential difference be¬ 


tween the up and down spin electrons in wire k. If the 
TLL wires are connected to Fermi-liquid leads, then the 
charge and spin conductance tensor is 

G£ ll = ^-(I-M„) . (12) 

If there are no Fermi liquid leads, and the TLL wires 
extend to infinity, then the charge and spin conductance 
tensor is given by 

Gtll ^ 2 ^® 2 (h-m„) . (13) 

We note here that Eq. (12) and Eq. (13), are strictly 
valid in equilibrium and at zero temperature. Any devi¬ 
ation from the equilibrium (say arising from a finite bias 
voltage or finite temperature) will lead to g dependent 
power law corrections, either on temperature or on the 
bias voltage or the system size, depending on whichever 
corresponds to the shortest length scale in the problem 
as in the case of a two wire junction [ 12 ]. 

Having described the various fixed points at the junc¬ 
tion and the associated conductance for both the Fermi 
liquid and TLL leads, we now proceed to calculate the 
TDOS of the spin-1/2 Y-junction. 


III. TUNNELING DENSITY OF STATES 


In this section we compute the TDOS of a spin-1/2 Y- 
junction for adding an electron with energy tux in one of 
the wires. The expression of the total TDOS in i th wire 
is given by the sum of the spin resolved TDOS Pi(tx) = 
Pfi(w) + A|,i(w), where 

i r°° / \ 

PaiM = — y )}^cft. (14) 


The spin resolved Green’s function in the i th wire is G = 
(V’ t Tz(^ 5 i)^ii( a; ; 0))- For a system with long wires, i.e., in 
the L —> oo the Green’s function can be written as G = 
G/ + G 0 = (i> I ai (x,t)ilj I J i (x,0)) + (tOM)Vb?0u 0 )), 
where we have neglected two oscillatory terms which van¬ 
ish in the L —» oo limit, and are unimportant. The non 
oscillatory terms in the spin resolved Green’s function 
are explicitly given by 


0)) = <^(M)V£(z,0)} 

n 


(15) 


1 

27TO 


—ia 

9, 

— a 2 — 4x 2 

v„t — ia 


(v v t — ia) 2 — 4x 2 


where we have defined f3 v = (1 + gD/i^g^), 7 ^ = (1 — 
gl)/(8gv) and dvi denotes the i th diagonal elements of 
the corresponding M„ matrix representing the boundary 
condition at the junction and v = c/s. For fixed points 
belonging to the M lc M ls class of fixed points d ui = a ui 
and for fixed points belonging to the M 2 C M 2 S class of 
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fixed points, = b u , c„, a„ for wires one, two, and 
three, respectively. Note that far away from the junction, 
in the x —> oo limit, the last two terms in Eq. (15) become 
unity and do not contribute to the TDOS, while in the 
x — > 0 limit, they contribute to the TDOS making the 
TDOS exponent different at the junction [41]. 

Inserting Eq. (15) in Eq. (14) and performing the in¬ 
tegration we obtain analytical expressions for the TDOS 
for the limiting cases of x —» 0 (in vicinity of the junc¬ 
tion) and x — > oo (far away from the junction). In both 
cases, the TDOS for w > 0 has the same analytical form 
and is given by 

*•<“> = . W 

where T(x) represents the gamma function, r c = a/v is 
an effective short time (inverse high frequency) cutoff, 
and uj = eV/h , where V is the bias voltage between the 
STM tip and the i th wire [see Fig. 1], In Eq. (16), A,; 
is the spin independent TDOS exponent which can be 
expressed as a combination of two terms: A* = A ci + 
A s ,;, where A c j (A s j) is only a function of g c ( g s ) and 9 C 
(9 S ) which specifies the boundary condition of the charge 
(spin) sector at the junction. 

Note that the TDOS exponent in general depends on 
whether it is being measured far away from the junc¬ 
tion, in the vicinity of the junction or at intermediate 
locations. Below we describe each of these regimes in 
subsections III A, IIIB, and IIIC, respectively. 


A. TDOS exponent far away from the junction 

Far away from the junction, i.e., in the x —> oo limit, 
for all possible classes of fixed points we find that 

A ” = i( ft + y) + i( 9 ' + y) ■ (17) 

independent of the boundary condition (9 C / S ) at the junc¬ 
tion, as expected. Equation (17) also specifies the TDOS 
exponent for a single infinite wire made of spin-1/2 TLL 
and is well known in the literature [41, 48, and 52]. We 
emphasize that even for the case of g s ^ 1 where the 
SU2 symmetry is broken in the wires, the TDOS expo¬ 
nent is identical for the up-spin and down-spin electron 
tunneling since the Hamiltonian for each wire in Eq. (6), 
and the boundary conditions at the junction are invari¬ 
ant under the interchange of spins, i.e., f-spin —> ./-spin 
and vice versa. Note that the result for the spineless case 
can be reproduced by substituting g c —)• g and g s —> g in 
Eq. (17). 


B. TDOS exponent in the vicinity of the junction 

In the vicinity of the junction, x —> 0, A, depends 
on the boundary conditions of both the charge ( 6 C ) and 



FIG. 2. Contour plot of the TDOS exponents in the vicinity 
of the junction, Ao, Ai, An and A 3 in the 6 C — g c plane for 
the SU(2) symmetric line g s = 1, in panels a), b), c), and d), 
respectively. Note that for g a = 1, the TDOS exponents in 
Eqs. (18) and (19) are independent of 9 S . The fixed points 
corresponding to the dome shaped regions lying to the left of 
the g c = 1 vertical line in all four panels show an anomalous 
enhancement of the TDOS, even for repulsive interactions. 
The boxed numbers specify the numerical value of the expo¬ 
nent for the corresponding contour lines. 


the spin ( 9 S ) sectors in general. For boundary conditions 
where both the charge and the spin sectors belong to the 
Mi c Mi s class, the Y-junction has Z3 symmetry in the 
wire indices and the TDOS exponent is identical in all the 
wires. The TDOS exponent is given by A 0 = A 0c + A 0s , 
where 

A 1 l + 5gg + (flg-l)cosfl„ 

°" 6 g v 1 +gl + {gl~ 1) cos 9 V 

where v = c/s. As a consistency check we note that 
as g v —> g , and 9 V —> 9 , the TDOS exponents for the 
spin-1/2 TLL Y-junction become identical to the case of 
spinless three wire junction, reported in Ref. [33]. 

For the case of broken Z 3 symmetry, i.e., where both or 
either of the charge and spin sector boundary conditions 
at the junction are specified by the M2 class of matrices, 
the TDOS exponent, in the vicinity of the junction, ex¬ 
plicitly depends on the wire index i. Particularly for the 
M2 C M2 s class of fixed points, the exponent is given by 
A i = A i c + A i s , where the spin and charge part of the 
exponent for wire 1 are explicitly given by 


_ 4 + 2s£ + (cos 9 V - vising/)(^ _ 1 ) 

— -ps- ’ 

12 g„ 

For the other two wires, A 2i , and are simply ob¬ 
tained by shifting 9 V -A 9 V q= 27 t/ 3 respectively, in the 
corresponding expressions for A 3 ^ in Eq. (19). For the 
Mi c M 2 S class of fixed points, the TDOS exponent is given 
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by Aj = Ao c + Aj S , and for the M 2 S Mi c class of fixed 
points, it is given by Aj = Aj C + Ao s , where Ao„ is given 
in Eq. (18), and Aj„ is specified by Eq. (19). 

For the particular case of the SU(2) symmetric line 
g s = 1, the spin part of the TDOS exponent in both 
Eqs. (18) and (19) become s A j s = 1/2, independent of 
9 S , where j = 0,1,2 or 3. We plot the TDOS exponent 
for this particular case, in the 9 C — g c plane in Fig. 2. 
In all four panels, the dome shaped region to the left of 
the g c = 1 vertical line, indicates fixed points showing 
an anomalous enhancement of the TDOS for small bias 
voltage even for repulsive interactions, in the vicinity of 
the spin-1/2 Y-junction. 

To gain more insight into the behaviour of the fixed 
points corresponding to the TDOS enhancement in the 
weekly interacting region ( g c ss 1) for the g s = 1 line, 
we expand the TDOS exponents to the lowest order in 
the small parameter (1 — g c ) and obtain, Aj = 1 + (1 — 
gc)dj C / 2, where dj C = a c , b c , c c or a c for j = 0,1,2 or 3, 
and it simply denotes the diagonal elements of the corre¬ 
sponding M c matrix, which relates the non-interacting 
incoming and outgoing charge fields at the junction. 
Thus, for weekly repulsive interactions, g c < 1, we find 
that Aj is less than one (or TDOS enhancement occurs) 
wherever dj C < 0, which physically corresponds to a hole 
current being reflected from the junction for an incident 
electron current. This is consistent with the previously 
reported enhancement of the TDOS for a TLL wire con¬ 
nected at one end to a superconductor [53] , where a prox¬ 
imity effect induced Andreev process leads to the reflec¬ 
tion of a hole for an incident electron. We emphasize here 
that even though we do not have any superconductivity 
explicitly in our model, the current splitting matrices, 
span all possible scenarios and includes the cases where 
a hole (either fully or partially) is reflected at the junc¬ 
tion. 

Let us now consider some specific fixed points for the 
three wire junction of spin-1/2 electrons, which have a 
unique basin of stability in the g c — g s plane [31], starting 
with the fixed points belonging to the Mi c Mi s class. For 
the NN fixed point Aq = 1/(2 g c ) + 1/ (2g s ) , and for g s = 1 
the range of g c where Ao < 1, or TDOS enhancement 
occurs, is given by g c > 1. Since the NN fixed point rep¬ 
resents three disconnected wires, this result is consistent 
with earlier results for both spin-1/2 [41] and the spinless 
disconnected wires [5, 49-51] where the TDOS enhance¬ 
ment at the ends of a single TLL wire can only occur for 
attractive interactions ( g c or g > 1). For the DD fixed 
point ( 9 U = 7r), A 0 = (1 + 2g 2 )/(6g c ) + (1 + 2g 2 )/(6g s ), 
which implies that for g s = 1 the TDOS enhancement 
occurs for 1/2 < g c < 1, i.e., even for repulsive inter¬ 
actions. For the DN fixed point ( 9 C = tv and 0 S = 0) 
A 0 = (1 + 2g 2 )/(6g c ) + 1/(2 g s ), and like the DD fixed 
point, the TDOS enhancement for the SU(2) symmetric 
line g s = 1 occurs for 1/2 < g c < 1. In case of the ND 
fixed point ( 9 C = 0 and 9 S = tv), the TDOS exponent is 
A 0 = 1/(2 g s ) + (1 + 2g 2 )/(6g s ), and for g s = 1 TDOS 
enhancement only occurs for attractive interactions as 






gc 


gc 




gc 


gc 


FIG. 3. Contour plot of the TDOS exponents in the g 3 — g c 
plane for various fixed points. Panels a), b), c), d), and e) 
display Ao for various fixed points of Mi c Mi s type, i.e., NN 
( 9 C = 9 S = 0), DD ( 9 C = 9 S = 7r), DN ( 9 C = tv,9 3 = 0), ND 
( 9 C = 0 ,9 S = 7r), and X + X + (9c = 9 3 = 2vr/3) respectively. 
Panel f) denotes TDOS exponent Ai for the DaDa fixed point 
of the M 2 C M 2 s type. The region bounded by Aj < 1, in 
all the panels signifies TDOS enhancement in the zero bias 
limit and we clearly see enhancement in the TDOS even for 
repulsive e-e interactions ( g c < 1) in several cases. The boxed 
numbers specify the numerical value of the exponent for the 
corresponding contour lines. 


in the NN case, i.e. for g c > 1. For the chiral fixed 
points, x + X + and X~X~, A 0 = (1 + 3g 2 c )/(2gl + 6 g c ) + 
(l + 3g 2 )/(2g 3 +6g s ), and the TDOS enhancement again 
occurs only for g c > 1 for g s = 1. 

Another fixed point with a unique basin of stability is 
the DaDa fixed point, which belongs to the M 2 C M 2 S class 
and is Z 3 asymmetric in the wire indices. Let us specif¬ 
ically consider the case of 9 C = 9 S = 0, in which wires 
1 and 2 are fully connected in both the charge and spin 
sectors and wire 3 is completely disconnected. For this 
case Ar = A 2 = Aoo = (1 + g 2 )/(4g c ) + (1 + g 2 s )/(±g s ), 
as expected in the bulk of an infinite TLL wire, and 
A 3 = 1/(2 g c ) + 1/(2 g s ) as expected, consistent with 
TDOS exponent near the ends of a TLL wire in Ref. [19 
and 41]. However, in this case, there is no anomalous 
TDOS enhancement either in wires 1 and 2 for any value 
of g c , or in wire 3 for repulsive interactions, along the 
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g s = 1 line. In Fig. 3, we plot the spin independent 
TDOS exponent in wire 1 (Ai) in the g s — g c plane 
for several fixed points (NN, DD, DN, ND, x + X + and 
D^D^) and clearly see regions where the TDOS is en¬ 
hanced (Ai < 1) in the zero bias limit, even for repulsive 
e-e interactions. 

Note that the D fixed point considered above also 
describes the two wire junction R^Rd where both charge 
and spin sectors are connected (wires 1 and 2), and the 
RnRn case where both the charge and spin sectors are 
disconnected (wire 3). The case of RdRn fixed point 
where the charge sector is connected and the spin sec¬ 
tor is disconnected can be constructed by considering 
wires 1 and 2 in the M 2 c Mi s class with 9 C = 9 S = 0 
(also called the D^N fixed point in the three wire con¬ 
text) and the TDOS exponent in this case is given by 
Ai = A 2 = (1 + 5 c)/( 4 5c) + 1/(2 g s ), which for g s < 1 
can never have TDOS enhancement for any value of g c . 
Finally, the case of RnRd fixed point with the spin sec¬ 
tor connected while the charge sector is disconnected, 
is equivalent to considering wires 1 and 2 in the ND^ 
fixed point (of Mi c M 2s class) with 9 C = 9 S = 0, and in 
this case Ai = A 2 = 1/(2 g c ) + (1 + <?s)/(4g s ), which for 
g s = 1 can have TDOS enhancement only with attractive 
interactions (g c > 1). 


C. TDOS crossover from boundary to bulk at 
finite distance from the junction 


At finite distance from the junction, x ^ 0, TDOS can 
be obtained from Eq. (14) after substituting Eq. (15), and 
performing the integration over t in the upper half com¬ 
plex plane, along the five branch cuts running parallel to 
the imaginary axis with branching points ia, ia±2x/v c , 
and ia ± 2x/v s . A similar approach is described in Ap¬ 
pendix A of Ref. [54] where the authors studied TDOS 
in a spiral TLL wire in proximity to a superconductor. 
For 0 < <jj < hv m i n /a , where u m i n = min{u c ,u s }, only 
the first term in Eq. (14) contributes, and the TDOS in¬ 
tegral in Eq. (14) becomes equivalent to the sum of five 
contour integrals around each of the branch cuts. The 
TDOS asymptote at finite x > v mm /(2ui) is obtained to 

be p ia (x,u) « Poo + pfa + Pia > where p x does not de¬ 
pend on the spatial coordinate and it is given by Eqs. (16) 
and (17). For x > v m i n /(2to), the other two terms dis¬ 
play an oscillatory power law behaviour on the spatial 
coordinate and are explicitly given by 


O) 

Pj 


(x,w) 


21-Aoo v Pb+i pdm 

'Kv v T{^ v d ui )v^ (y v - v D )^dui 



x 



/ 2 xuj 
V 



( 20 ) 


where 8 V = Arg(z Aoo+7,/<i ‘' i ), (u, v) = (c, s) or (s, c), and 
/3„, 7 „ and d j„ are defined below Eq. (15). Along the 
SU(2) symmetric line, g s = 1, = 0 and the depen¬ 

dent term in Eq. (16) drops out of the TDOS integral, 


and the spatial dependence of the TDOS arises only from 
the charge term: pf 1 oc cos(2xu}/v c — 5 c )x 7c ^ oi_/3c_1 / 2 
with a slowly oscillating contribution that drops off as 
a power law. Note that in the vicinity of the junction 
[x < u m i n /(2w)], the TDOS has a power law behaviour: 
Pi oc w Ai_1 where A* depend on the boundary con¬ 
ditions at the junction [see Eqs. (18) and (19)] . For 
small distances away from the junction [a; ~ u m i n /(2w)], 
TDOS shows an oscillatory behaviour with bias fre¬ 
quency and for large distances away from the junction 
[x iViin/(2a>)], the power law behaviour in x [see 

Eq. (20)] suppresses the oscillations and the TDOS goes 
over to the bulk value, pi oc uj A oo ~ 1 . 


IV. CONCLUSION 

To summarize, in this article we explicitly calculate 
the local tunneling density of states in the vicinity of 
a spin-1/2 TLL Y-junction, and present an analytic ex¬ 
pression for the TDOS exponent in terms of the boundary 
condition at the junction and the strength of e-e inter¬ 
action. We find that there are several fixed points which 
in the vicinity of the junction give an anomalous TDOS 
enhancement in the zero bias limit, even for repulsive 
interactions. Physically all such instances of TDOS en¬ 
hancing fixed points are such that holes are reflected from 
the junction for incident electrons. This makes the TLL 
junction physically similar to the case of a TLL connected 
to a superconductor where the TDOS enhancement was 
attributed to the proximity induced Andreev processes 
[53], even though the three wire junction considered by 
us has no superconductor. 

It should be noted that the TDOS expression in the 
vicinity of the junction, given by Eq. (16), is valid only 
for x < x c = v m i n /(2ui), where u m i n = min{u c , Vo-}. For 
a realistic system, u m i n ~ 10 5 m/s, and for a STM tip 
voltage of lmV, the crossover length scale is x c ~ 33 
nm, which is within current experimental reach. For x ss 
x c , the TDOS displays an oscillatory behaviour which 
is again suppressed at large distances, x x c , and the 
TDOS reverts back to its bulk value. One limitation of 
our calculations is that they are valid only for the regime 
of bias frequencies which do not breach the linearization 
regime of each TLL wire, i.e., w < u m i n /a:. 

Note that while we have considered a junction with 
symmetric e-e interaction strength, in a more realistic 
situation this may not be the case due to asymmetrical 
screening induced by gates, inhomogeneities, or defects. 
A detailed analysis of TDOS in a Y-junction with differ¬ 
ent values of g in the three wires, can be done in the spirit 
of Ref. [38], which focused on the stability and analysis 
of the fixed points. However it is beyond the scope of the 
present paper, and can be the subject of a future work. 
Additionally, we note that all our results are valid only 
in the regime where backward and umklapp interactions 
can be safely ignored [55]. 

In addition to the spin degree of freedom, our TDOS 



exponent calculations can also be extended to include 
other degrees of freedom such as different valleys in 
carbon nanotubes. If the spin and the valley rotation 
symmetry are not broken, then all the corresponding 
TDOS exponents for carbon nanotube Y-junctions are 
easily calculated using the fact that three of the four 
bosonic fields get pinned and only the center of mass 
field which corresponds to the charge degree of freedom 
primarily contributes interaction dependent term in the 
TDOS. The TDOS exponents for carbon nanotube Y- 
junctions are explicitly given by A^ N = (2A ]C + 3)/4, 
where j = 0,1,2 or 3, and A JC are specified by Eqs. (18)- 
(19). This is consistent with the TDOS exponents re¬ 
ported for the bulk and the edge of a TLL hosted in a 
carbon nanotube [55]. 

Experimentally, spin-1/2 TLL wire junctions can be 
realized by means of carefully patterned ID wires in a 
2DEG, or via crossed single-walled nanotubes [21] and 
tuned to various fixed points by means of nano-gates ap¬ 
plied near the junction. Another feasible possibility is an 
island like set-up of quantum Hall edge states, proposed 
in Ref. [28] , in which the tunneling operators can be con¬ 
trolled via gate voltage operated constrictions in the cen¬ 
tral region of the island. Once the junction is tuned to 
an appropriate fixed point, the power law exponent of 
the TDOS can be experimentally extracted by measur¬ 
ing the differential tunneling conductance as a function 
of the STM tip voltage (for fixed temperature), or as a 
function of temperature (for fixed voltage) as was done 
in Ref. [19]. 
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Appendix A: Scaling dimensions of the fixed points 
at the junction 

In this appendix, we present the scaling dimensions of 
all the spin-preserving single particle and pair tunneling 
operators, for all possible fixed points. For any boundary 
operator O b , the scaling dimension Sq can be calculated 
by the two point correlation function 

, (Al) 

and it depends on the strength of the e-e interactions 
(< 7 „), and the boundary condition at at the junction ((9„). 


The stable boundary condition or fixed points are those 
for which the scaling dimension of all boundary opera¬ 
tors is either <5o = 0 or S 0 > 1 for a given e-e interaction 

TABLE I. Scaling dimension of various single particle tunnel¬ 
ing operators: <5oo- = 8kc + Si a , where fc, l = 1 or 2 and the 
boundary condition at the junction is specified by Mfe c Mi s . 


Operator 

Slv 


2g u (1 — cos 9 U ) 
3[g^+cos0„(g^-l) + l] 


g u ( cos 9 l/ -\-y/3 sin 9 u -\-2) 

3[g;?+cos0„(g2_l) + l] 

V’gaV’l.a, 

g u (cos 9 U — \/3 sin 9 U +2) 
3[g;f+cos0„(gj; —1) + 1] 


9u 

gj^+cos 9 U (g‘£ — 1)+1 

V'ga» V4V ^2,0-, *P°<T 

9v 

g^ -|- cos 9 U (g^ — 1) + 1 

Operator 

$2i/ 

i’lji’lcy 

Ig„(2 - 2cos6b) 

i>°,J ipl* 

!<?„( 2 + cos#„ + \/3sin6b) 


|<?„( 2 + cos 6 V — \/3sin6b) 


|^(2 -2 cos 6b) 

V’3,<r,'/'°aV’2, CT 

^^j(2 + cos 6b — \Z3sin0„) 


H^(2 + cos 6b + \/3sin6b) 

2,<r 

(s2+l)+2(aS-l)coBfl„ 

4 9„ 


(< 7 ^ + l) — (g„ — l)(cos 9 U — \/3 sin 9 U ) 

4 s„ 


(g^+1) — (g^ — l)(cos 9 v -\-y/ 3 sin 9 U ) 

4 9u 


strength. For a fixed point specified by Mfc c M; s , where 
k,l = 1 or 2, the scaling dimension of all single particle 
and pair tunneling operators can be expressed as a sum 
of the charge and spin components: <5o = 5k c + $is- For 
the single particle tunneling operators, 6 i„ is explicitly 
given in the upper half of Table I, and is tabulated 
in the lower half of Table I. The scaling dimensions of 
all possible pair tunneling operators, is tabulated in Ta¬ 
ble II. Finally, we note that, ideally we need to know the 
scaling dimensions of all possible multi particle tunneling 
operators to determine the stability of a fixed point, how¬ 
ever, more particle processes tend to be less relevant, and 
based on a conformal field theory argument it has been 
explicitly shown in Ref. [31] that single and two particle 
tunneling operators completely determine the stability of 
all the fixed points which have a unique basin of attrac¬ 
tion for a spin-1/2 Y-junction. 
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TABLE II. Scaling dimension of various pair tunneling operators: Soc = Skc + Si s , where k,l = 1 or 2 and the boundary 
condition at the junction is specified by the matrix Mfe c M; s . 


Operator 

Sic 

<5i s 

$2 c 

<52s 


4g c (2+cos 0 c + \/3 sin 0 C ) 

0 

2 i|Jfi-cos 0 c ) 

0 

3(l+g2_|_( fl ,2_ 1 ) cos 0 c ) 


?) 

0 

(2 + cos 9 C — VS sin 0 C ) 

0 



0 

3 g g 3 ° (2 + cos 0 C + sin 9 C ) 

0 


49c (2+cos 0 C — \/3 sin 0 C ) 

0 

+ cos ^c + VS sin 0 C ) 

0 

3(l+92 + (g, 2 _i) cos 9 c ) 


?) 

0 

^(i-cos 0 c ) 

0 



0 

(2 + cos 0 C — VS sin 9 C ) 

0 

V’ 2 ,aV’l.aV'° CT + V’°a 

29 c(l+c°s 0 C ) 

29 s ( 1 +cos 0 S ) 

9(1—cos 0 C ) — \/3 sin 0 c (l+2 cos 9 C ) 

9(1 —cos 0 S ) — \/3 sin 9 S (1+2 cos 0 S ) 

l+9? + (9?-l)co S e c 

l+sf + tsl-l) cos 9 S 

9 9c 

9 9s 

V’3,JV’laV'° CT + V’°a 


?? 

9(2+cos 0 C ) — \/3 sin 0 C (11+4 cos 0 C ) 

9(2+cos 0 S ) —\/3 sin 0 S (11+4 cos 0 S ) 



189c 

189s 



V 

9(2+cos 0 c ) + y/3 sin 9 C {7 —4 cos 0 C ) 

9(2+cos 9 S ) + >/3 sin 0 S (7—4 cos 0 S ) 



18 9c 

189, 


29 c( 1 +cos 0 C ) 

29 S ( 1 —cos 0 S ) 

9(1—cos 9 C ) — V3 sin 9 C (1+2 cos 9 C ) 

9 S (9 (1+cos 9 S ) — \/3 sin 9 S (1+2 cos 0 S )) 

l+9c+(9c- 1 ) cos0 c 

l+s»+(gj}-l) cosSs 

9 9c 

9 

v , 3 ,iV' 2 , CT v , °- t CT V'£- CT 


5? 

9(2+cos 0 C ) — \/3 sin 0 c (ll+4 cos 0 C ) 

9 S (9(2—cos 0 S ) + \/3 sin 9 S (7—4 cos 0 S )) 



189c 

18 

V’UV’laV'p-aV’^-cT 


» 

9(2+cos 0 c ) + \/3 sin 9 C ( 7—4 cos 0 C ) 

9 s (9(2—cos 9 S ) — \/3 sin 0 S (11+4 cos 0 S )) 



189c 

18 


89 c ( 1 —cos 0 C ) 

0 

+p (2 + cos 9 C — \/3 sin 9 C ) 

0 

3(l+9c + (fl'c —1 ) cos 0 c) 

^24^14^24^24 

?) 

0 

(2 + cos 0c + %/3 sin 9 C ) 

0 

V'°t t V’34V , °/V’34 

?) 

0 

%+l - cos 9 C ) 

0 

,-<7 

29 c ( 1 —cos 0 C ) 

29 S (1 —cos 0 S ) 

^(1 — cos 9 C ) 

g s ( 1 + cos 0 S ) 

3(l+Se+(9c-l) cos ®c) 

(l + <)2 + (g 2 _l) cos9s ) 

^2^ ^2,^3,-^3 ,- CT 



(2 + cos 0 C — v^3 sin 0 C ) 

^ (2 — cos 0 S + VS sin 0 S ) 

V'° CT t V’3,a-*/’u- t CT V ) l,- CT 


V 

2 a (2 + cos 0 C + \/3 sin 0 c) 

^ (2 — cos 0 s — VS sin 0 3 ) 


29 c ( 1 —cos 0 C ) 

29 S ( 1 —cos 0 S ) 

^(1 - cos 0 C ) 

^(1 - cos 0 S ) 

3(l+9§ + (g- 2 _i) cos 9 c ) 

3(l+g| + (g2 — 1 ) cos 0 S ) 

i>°J V>3,<t 

?) 

» 

^ (2 + cos 0 c — VS sin 0 C ) 

^ (2 + cos 0s — \/3 sin 0 S ) 

V’° CT t V'3, CT V’n CT Vl,a 



(2 + cos 0 C + \/3 sin 0 C ) 

^ (2 + cos 0 s + \/3 sin 0 3 ) 

V’°a t V’l,aV’°- t CT V’ 2 ,- CT 

29 c ( 1 —cos 0 C ) 

29 S ( 1 +cos 0 S ) 

(1 — cos 0 C ) 

^(1 - COS 0s) 

3(l+9§ + (g- 2 _i) cos 0 c ) 

(l+ 9 f + ( 9 l —l)cos 0 s ) 

V’° CT t V’2,aV , °- t CT V’3 ,-CT 

» 

5? 

^ (2 + cos 0 C — \/3 sin 0 C ) 

2 ^— (2 + cos 0 s — v^3 sin 0 3 ) 

V'p CT t V’ 

?) 

» 

(2 + cos 0 C + VS sin 0 C ) 

2 ^ 7(2 + cos 0 s+ ^3 sin 0 3 ) 

V’°t t ^btV'°tV'[ 4 t 

0 

49 s (2+cos 9 S + \/3 sin 0 S ) 

0 

^( 1 -cos 0 3 ) 

3(1+9^ + (9^-1)cos0„) 

V'°t t ^2.t 1 / , 3!t^24 t 

0 

5? 

0 

;i ^ i (2 + cos 0 s - VS sin 0 3 ) 

V’°t t ^34V’b4.V'34 t 

0 

» 

0 

(2 + cos 0 s + \/3sin0 s ) 

V’°t t ^btV'°tV'[4 t 

0 

49 s (2+cos 9 S — -\/3 sin 9 s ) 

0 

/ O I p/~»c ft I « /q C 1 T 1 ft 'i 

3(l+g| + ( 9 g—1) cos 0 S ) 

Qg s V" • COb Us 1 V bill (7s ^ 

V’ut^24' ! /’u4-^ , 24 

0 

V 

0 

2 ^( 1 - cos 0 s) 

V’°t t ^34V , °tV , 34 t 

0 

» 

0 

(2 + cos 0 C — VS sin 0 3 ) 

V’°a t V'l.<7'0°-aV , 2 

29 c( 1 +cos 0 C ) 

29 s (1—cos 9 S ) 

^(1 - cos 0 C ) 

^■(1 — cos 0 3 ) 

1 + 92 + ( 9 2 - 1 )cos 0 c 

3(1+9, +(9s“l) cos 6 S ) 

V’s^V’i 2,+°-<+3 ,-<t 

?) 

V 

2 ^- (2 + cos 0 C — VS sin 0 C ) 

(2 + cos 0s — \/3 sin 0 3 ) 



V 

2 ^ (2 + cos 0 c + VS sin 0 C ) 

7 + (2 + cos 0s + \/3 sin 0 3 ) 
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